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ABSTRACT 


This report presents the results of the first phase of a three-phase study on the 
development and application of the energy-based method for prediction of the 
liquefaction potential of sandy soils. The formulation of the method is based on the 
convolution of the basic elements from both the “stress” and “strain” approaches and is 
very flexible in incorporating the special characteristics of ground motion such as the 
near-field effects. The feasibility phase consists of the tasks: 1) to collect and synthesize 
laboratory data; 2) to perform ground response analyses at the Wildlife Site, which 
suffered a massive ground liquefaction failure during the Superstition Hills Earthquake; 
and finally 3) to compare and to assess the differences between the field and the 
laboratory data. Even though the scope of the feasibility study did not permit cyclic 
testing of the soil samples from the Wildlife Site, the correlation of the field response 
data and the applicable laboratory data are strong. The results of this phase suggest that 
development of an energy-based method to evaluate liquefaction potential is feasible. 


KEYWORDS: building technology; liquefaction; strain energy; earthquake; ground 
response; cyclic testing; laboratory measurements; ground motion; pore pressure. 
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CHAPTER 1 
INTRODUCTION 


11 BACKGROUND 


Liquefaction failure has been and continues to be a major cause of damage during 
earthquakes. The direct and indirect costs associated with ground failure may far exceed 
the damage caused by other types of failures such as structural collapses. Due to the 
enormous damage potential, research in the areas of liquefaction prediction and 
mitigation has continued, and the respective technologies have significantly improved 
over the years. 


Two basic methods are currently used to predict liquefaction potential. The most widely 
used method, based on laboratory data and field performance data, was developed by 
Seed et al. (1983, 1985). In this method, the cyclic stress ratio in the field is predicted 
based on a simple or a more detailed ground response analysis, and the demand resulting 
from the design earthquake is established. The cyclic shear strength of the material can 
be obtained from laboratory testing of the soil samples or from the penetration data 
(Standard Penetration Test [SPT] data or Cone Penetration Test [CPT] data) along with 
the index and gradation properties of the soil samples. The most widely used curves to 
predict the capacity of the soil in terms of cyclic shear strength are the set of curves by 
Seed et al. (1985), shown in Figure 1.1. From knowledge of the penetration resistance 
and the fines content in the soil materials, the cyclic shear strength can be determined. 
However, the application of this method involves several empirical factors, including the 
corrections for sample disturbance, earthquake magnitude, and overburden pressure. 
Recently, Arango (1994) presented the new developments for this method, including the 
effects of higher frequency of loading on the cyclic shear strength of the soil and most 
notably, the recommendations made for revising the correction factors for earthquake 
magnitudes. Another recent study (Koester, 1992) has shown that the correction for 
overburden pressure is also a function of the fines content and that a reduction of the 
cyclic shear strength due to overburden pressure will significantly decrease if the fines 
content of the materials increases. The recent publication by Ishihara (1993) also 
provides an adjustment factor to incorporate the plasticity effects of the fines on the 
cyclic shear strength of the materials. Over the years, the “stress” method has proved to 
be a conservative and reliable method for the prediction of liquefaction potential, 
especially for distant earthquake events, which was the basis of the data used for 
development of this method. However, the method lacks the flexibility to incorporate 
recently recognized characteristics of the earthquake ground motions such as the near- 
field effects. For example, the “fling” effect resulting from the source and directivity of 
the rupture which was recently observed in the Kobe and the Northridge earthquakes, 
concentrates most of the energy in a short period of time. Such effects, combined with a 
much higher intensity of the ground motion in excess of 1g peak ground acceleration 


recorded in urban areas in recent major earthquakes, require a more robust approach to 
investigate the liquefaction potential in such regions. 


The second method for predicting liquefaction potential is based on the “strain” 
approach. This method was developed among others by Dobry et al. (1982). In this 
method, the shear strain in the field is compared with the laboratory data relating cyclic 
shear strain to excess pore pressure to determine the liquefaction potential. Similar to the 
“stress” method, the “strain” method also requires ground response analyses and 
laboratory testing of the soil samples. The “strain” method is fundamentally different 
from the “stress” method and lacks the wide range of the field and the laboratory data 
bases that exist for the “stress” approach. 


The “strain energy” method discussed in this report incorporates the basic elements of 
both the stress and strain approaches in the formulation. In this method, the amount of 
total strain energy at the onset of liquefaction is obtained from the stress and strain time 
histories from laboratory testing and is compared with the same energy in the field due to 
the design earthquake motion. The basis for this method is the observation made on the 
laboratory data that the build-up of the excess pore pressure is proportional to the total 
strain energy in all loading cycles up to the initial liquefaction. This observation has 
prompted the formulation of the “energy-based” method. This method has been 
investigated in recent years by several researchers, including Figueroa et al. (1994, 1995) 
and Kagawa et al. (1990). 


1.2 PURPOSE 


The purpose of this study is to evaluate the feasibility of the development and application 
of the strain energy method for general use. The study is expected to continue with two 
additional phases that will develop generic “strain energy” liquefaction curves as a 
function of the most relevant soil properties and generic “strain energy” demand as a 
function of seismicity data and a wide range of site soil data and profiles. The limited 
scope of the feasibility study did not permit laboratory testing for the purpose of the 
“strain energy’ computation. Available laboratory data were used for this purpose. 


13 OVERVIEW OF THE REPORT 


In this report, Chapter 1 includes the introduction and scope of the study. Chapter 2 
presents the collection and synthesis of the laboratory data. Chapter 3 discusses the soil 
and earthquake data from the Wildlife Site. The ground response analyses and 
comparison of the results with the laboratory data are presented in Chapter 4. Finally, 
Chapter 5 presents the summary and the recommendation. The references are listed in 
Chapter 6. All the laboratory data used in this report are presented in Appendices A 
through D. 
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Figure 1.1 - Relationships Between Stress Ratio Causing Liquefaction and (N;)60 Values 
for Silty Sands for M = 7-1/2 Earthquakes (Seed et al., 1985) 
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CHAPTER 2 


COLLECTION AND SYNTHESIS OF LABORATORY DATA 


Computation of the strain energy requires access to the stress and strain time histories 
from cyclic (triaxial or simple shear) tests in the laboratory. Such data are usually 
computer storage-intensive and are not maintained for a long period. However, an 
attempt was made to collect the available and reliable data to characterize the strain 
energy. Most of the data were obtained in connection with various recent Bechtel 
projects. The laboratory data used in this study are from: 


e The cyclic stress- and strain-controlled tests on Monterey No. 0 sand, performed at 
the University of California, Berkeley. 


e The stress-controlled tests on soil samples from the Savannah River Site (SRS), 
performed at the University of California, Berkeley. 


e The cyclic torsional shear tests on clean and silty sands, performed at the University 
of Colorado. 


e The cyclic triaxial tests on clean sands, performed at Wayne State University. 


e The summary of the laboratory data reduced to a set of relationships to compute strain 
energy, as developed by Figueroa et al. 


Altogether, a total of 150 cyclic test data sets have been processed. A limited number of 
these were excluded in the process due to peculiar stress and strain patterns and 
incompleteness of the respective time histories. The computation of the strain energy 
from each data set and a discussion on the validity of each group of tests follow. 


2.1 STRAIN ENERGY COMPUTATION 


In a typical cyclic laboratory test, the stress, strain and pore pressure time histories are 
recorded. Typical recorded time histories for a strain-controlled cyclic triaxial test are 
shown in Figure 2.1. Hysteresis loops can be developed from the shear stress and strain 
time histories. The hysteresis loops corresponding to the stress and strain time histories 
shown in Figure 2.] are shown in Figure 2.2. From the shear stress, T(t), and the shear 
strain, y(t) at time t, the time history of the total strain energy up to time t, E(t), is 
computed from: 


2-1 


T | 
E(T) = 3, Z (t).Ay (t) (2.1) 


where t is the time, >) is summation over the time increment At up to time t, T (t) is the 
average shear stress from time t to t + At, and Ay(t) is the shear strain increment from time 
ttot+At. The strain energy for each cycle of loading amounts to the area inside the. 
hysteresis loop. The computation of the instantaneous energy and its summation over 
time intervals were performed until the onset of the liquefaction, at which time the pore 
pressure ratio reached a value of unity. The summation of the energy at this time, Eig, 
was used as the measure of the capacity of the soil sample against initial liquefaction 
occurrence in terms of the strain energy. 


2.2 CYCLIC TRIAXIAL TESTS ON MONTEREY NO. 0 SAND, 
PERFORMED AT THE UNIVERSITY OF CALIFORNIA, BERKELEY 
(UCB) 


The data were prepared as part of the Bechtel in-house technical research led by Arango 
(1994). The gradation curve for Monterey No. 0 sand is shown in Figure 2.3. The tests 
were both stress- and strain-controlled. The samples were prepared at relative densities 
ranging from 40% to 60%, and the loads were applied at frequencies ranging from 0.10 
Hz to 20 Hz. All tests were conducted at a confining pressure of 100 kPa. More detailed 
information about the testing program and the testing apparatus may be obtained from 
Riemer (Riemer et al., 1994). A total of 20 tests from this group were incorporated in this 
study. A summary of the test data, including the computed total strain energy to the onset 
of liquefaction for each test, is presented in Table 2.1. The recorded stress, strain, pore 
pressure (in terms of the pore pressure ratio r,), and the computed time history of the 
Strain energy for each test are shown in Appendix A. In addition to the strain energy time 
history, the energy time history normalized to the total energy at the time of r, = 1, (Exiq), 
is also plotted and compared with the r, time history, e. g. see Page A-3. As shown in 
these plots, the normalized strain energy increase follows the pattern of the pore pressure 
ratio increase and, on the average, shows a very good agreement for all the tests at a wide 
range of frequencies and at all the relative densities tested. The agreement holds whether 
the data are obtained from the stress- or the strain-controlled tests. As stated earlier, this 
observation was the basis for formulation of the strain energy method. 


A summary of the results in terms of the total energy as a function of relative density is 
shown in Figure 2.4. As expected, the total energy to the onset of liquefaction increases 
as the relative density of the sample increases. It can also be observed in this figure that 
the scatter in the strain-controlled test data is less severe than the scatter in the data from 
the stress-controlled test results. 
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The strain energy for each test as a function of the frequency of loading is plotted in 
Figure 2.5, which shows a decreasing total energy as the frequency of loading increases. 
In this figure, the frequency of loading has a more pronounced effect on the total energy 
obtained from the stress-controlled tests than the strain-controlled tests. In addition, the 
strain-controlled tests require lower total energy to develop initial liquefaction as 
compared to the stress-controlled tests. It should also be noted that for a typical strain- 
controlled test, the pore pressure build-up takes place at a much faster rate in the first 
several cycles of loading. On the other hand, in the stress-controlled test, the rate of pore 
pressure build-up increases towards the end of loading cycles. This observation can also 
be made from the shape and size of the respective hysteresis loops. In the strain- 
controlled tests, the largest loops are the earlier loops, and they decrease in size as the 
sample degrades due to the pore pressure build-up. The opposite trend takes place in a 
stress-controlled test, as shown in Appendix A, e. g. Pages A-4 and A-16. 


2.3 CYCLIC TRIAXIAL TESTS ON SOIL SAMPLES FROM THE 
SAVANNAH RIVER SITE, PERFORMED AT THE UNIVERSITY OF 
CALIFORNIA, BERKELEY (UCB) 


The laboratory program for this group of tests was developed as part of one of the 
Bechtel projects for the Department of Energy (DOE) at the Savannah River Site (SRS). 
Subsurface conditions for the site under consideration are shown in Figure 2.6. The soil 
layers of primary interest were the Tobacco Road (TR3 and TR4) and the Santee 
formations. 


A comprehensive site investigation program was conducted at the site. Relevant average 
soil properties of each soil layer at the SRS site, shown in Figure 2.6 are summarized in 
Table 2.2. 


Most of the cyclic load tests were conducted on undisturbed soil samples from the 
Tobacco Road formation from depths of 16 m to 23 m. As shown in Table 2.2, this 
material has an average fines content of 23%, including 9% clay content (minus 2 micron 
particle size) and an average plasticity index of 25%. A typical gradation curve for the 
Tobacco Road Materials is shown in Figure 2.7. A total of 22 cyclic stress-controlled 
tests at 1 Hz were performed (Riemer and Seed, 1994). The confining pressure ranged 
from 200 kPa to 750 kPa. A summary of the test data and of the total strain energy for 
each test in this group is presented in Table 2.3. Notable characteristic of this group of 
tests is the large confining pressure used in the tests and the relatively large fines content 
in the soil samples tested. The plots of shear stress, strain, total energy and normalized 
energy, pore pressure ratios, and the hysteresis loops for this group are shown in 
Appendix B. As shown in this appendix, the increase of the normalized energy in general 
follows the pore pressure ratio increase up to the pore pressure ratio of one. A summary 
of the total energy as a function of the confining pressure is shown in Figure 2.8. As 
expected, the total strain energy is greater for the samples tested at higher confining 
pressures. For the same confining pressure, tests on samples having a higher dry density 
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resulted in the development of a larger total energy. This trend is similar to the trend 
observed in Figure 2.4 with respect to the relative density of the samples. 


2.4 CYCLIC TORSIONAL TESTS ON SOIL SAMPLES, PERFORMED AT 
THE UNIVERSITY OF COLORADO (UOC) 


The time histories for this group of tests were provided by Koester (1992 ). The test data 
were developed as part of the research work for a doctoral dissertation at the University 
of Colorado (UOC). Only the time histories from nine tests were available. The tests 
were performed using the stress-controlled hollow torsional simple shear test apparatus. 
The cyclic loading was applied at a frequency of 0.1 Hz. Both clean sands and sands 
with fines content up to 45% were tested. 


The silty sand samples were prepared with a density such that the void ratio of the sample 
matched that of the parent clean sand at the selected relative densities. The confining 
pressure in the tests ranged from 200 kPa to 300 kPa. A summary of the soil data and the 
test results in terms of the total energy is shown in Table 2.4, whereas the gradation curve 
is Shown in Figure 2.9. A more detailed description of the sample preparation and testing 
program can be obtained from Koester (1992). 


The time histories of the stress, strain, pore pressure ratio, total energy, and hysteresis 
loops for this group are presented in Appendix C. In general, the hysteresis loops in this 
group of tests start with a few narrow loops followed by one or two large loops before 
reaching the initial liquefaction stage, suggesting a sudden contraction and collapse of the 
samples. This behavior may have been the cause of the relatively low densities of the 
samples. The test results in terms of the total strain energy as a function of relative 
density for the clean sand are shown in Figure 2.10. The results show a relatively large 
scatter in the energy at low relative densities. The results of the silty sand samples as a 
function of the confining pressure are shown in Figure 2.11. These results show 
relatively less scatter in the data. 


2.55 CYCLIC TRIAXIAL TESTS ON SOIL SAMPLES FROM THE 
NORTHRIDGE SITE, PERFORMED AT THE UNIVERSITY OF 
CALIFORNIA, BERKELEY (UCB) 


The samples for this group of tests were prepared as part of the National Science 
Foundation (NSF)/Bechtel research work led by Arango (Arango and Migues, 1996). As 
part of the test program, a total of 8 reconstituted clean sand samples were prepared and 
tested in a stress-controlled cyclic triaxial test device. The samples were prepared at 
relative densities ranging from 35% to 90%. The gradation curves for two soil samples 
are shown in Figure 2.12. A summary of the test data and of the computed total strain 
energy is presented in Table 2.5. Time history plots are included in Appendix D. Test 
results in terms of the total strain energy as a function of relative density is presented in 
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Figure 2.13. As shown previously, the total strain energy increases as the relative density 
increases. 


2.6 CYCLIC TRIAXIAL TESTS ON CLEAN SANDS, PERFORMED AT 
WAYNE STATE UNIVERSITY (WSU) 


The summary results of 91 cyclic triaxial stress-controlled tests on clean sands was 
presented in the doctoral dissertation by Al-Khatib (1994). The tests were performed at 
Wayne State University (WSU). All tests were performed on clean sands consisting of 
Monterey No. 0 and Kasumigaura sand (K-sand). The gradation curves for the two sands 
are shown in Figure 2.14. The breakdown of the tests is as follows: 


e 28 tests on K-sand at low frequency with cyclic reversal loading (two-way cyclic 
loading) 

e 28 tests on Monterey No. 0 sand at low frequency with cyclic reversal loading 

e 25 tests on Monterey No. 0, low frequency and one-way loading 

e 10 tests on Monterey No. 0 with earthquake simulated loading using the El Centro 
and Taft records 


Time histories of the test data were not available; however, the total energy in terms of 
axial stress/strain has been reported by Al-Khatib (1994). The total energy was converted 
to the total energy in terms of the shear strain and shear stress attributes and are 
summarized in Table 2.6. The results in terms of the total strain energy as a function of 
relative density for the cyclic reversal loading cases are shown in Figure 2.15 which 
shows a similar trend to the one observed in the UCB data (see Figure 2.4). As shown in 
this figure, both the K-sand and the Monterey No. 0 sand have similar capacity in terms 
of total energy and consistently show an increase of the total energy with an increase in 
the relative density. The scatter in the data appear to be minimal. 


The results in terms of the cyclic one-way and two-way loadings are compared in Figure 
2.16. As shown, the two-way loading results in lower total energy capacity as compared 
to the one-way loading. This trend is consistent with the intuitive indication that soil 
resistance to liquefaction will be higher due to the less damaging effects of the one-way 
loading. Finally, the results of the two-way loading are compared with the earthquake 
loading in Figure 2.17. The earthquake loading results in the lower total energy. 
Altogether, the results of this group of tests appear to be more uniform with little scatter 
in terms of the total energy. 


2.7 SUMMARY DATA BY FIGUEROA et al. 
A series of torsional shear hollow cylinder tests were performed on both clean sand and 


silty sand by Figueroa et al. (1994, 1995). Samples from the Reid Bedford sand (clean 
sand) were tested at relative densities ranging from 50% to 70%. The silty sand from the 
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Lower San Fernando Dam (LSFD) were tested at relative densities of 57% to 92%. The 
gradation curves for both materials are shown in Figure 2.18. Each sample was 
successively tested at confining pressures of 41.4 kPa, 82.7 kPa, and 124.1 kPa. Actual 
data points for this group of tests are not available. However, the authors performed 
regression analyses of the test results in terms of total strain energy and identified the 
most relevant parameters affecting the results of clean sand and silty sand. Based on the 
test results, the authors recommended the following relationships (Figueroa et al., 1995): 


Clean sand 
Log dE = 2.062 + 0.0039 oe + 0.0124 D, (2.2) 
Silty sand 


Log dE = 2.529 + 0.00474 oo (2.3) 


where 6E is the total strain energy in Joules/m’, 6, is the effective confining pressure in 


kPa, and D, is the relative density in percent. The relationship for clean sand shows the 
confining pressure as one of the variables. However, the importance of this parameter is 
very small due to the small coefficient associated with this parameter in Equation 2.2. 


2.8 SUMMARY OF ALL LABORATORY DATA 


Based on the results of the five groups of tests outlined above, summary plots have been 
prepared to evaluate consistency between the various test groups. 


For clean sand, the results of tests on Monterey No. O performed at the UCB (20 tests), 
the data on clean sands from the WSU (81 tests), the data from the Northridge samples 
(8 tests) also tested at the UCB, and the data from the UOC (4 tests) are compared with 
the relationship by Figueroa et al. in Figure 2.19. As shown in this figure, except for the 
data from the UOC, the remaining groups show a quite consistent pattern of the rate of 
energy dissipation and of the total energy absorbed. The confining pressure used in the 
tests at the UOC was at least 2 to 3 times larger than the pressure used for the rest of the 
tests. Also, the differences in the shape and size of the sand particles may have 
contributed to some of the differences in the results. For relative densities in the range of 
40% to 70%, the data from UCB, WSU, and Figueroa et al. are in relatively good 
agreement. 


For silty sands, the results from the Savannah River Site (SRS) are compared with the - 
data from the UOC and the relationship by Figueroa et al. in Figure 2.20. The fines 
content in each group are: 28% for the samples from the lower San Fernando Dam 
(Figuero et al., 1995), 20% to 45% for soil samples tested at UOC (Koester, 1992), and 
the average 23% for samples taken from the SRS site. The Plasticity Index of the 


2-6 


materials in the groups also varies from 10 to 25%. Unfortunately, the confining 
pressures used for each group of tests do not overlap. Nevertheless, each group of results 
follows the pattern of the previous group and a consistent trend is maintained. 


All of the results indicate that for clean sands, the energy to liquefaction can be quantified 
in terms of the relative density and the confining pressure. However, the limited data 
available does not permit a study of effects of the grain size and shape on the total 

energy. 


Summary of the results for silty sands also shows that the energy to liquefaction can be 
quantified in terms of the effective confining pressure. However, the effects of the 
plasticity index, the amount, and the type of fine need to be studied in the future. 


As stated earlier in the report, the scope of this feasibility study did not include laboratory 
testing. However, comparison of the data available from the various researchers and 
practitioners at different institutes shows remarkably good agreement. This observation 
leads to the conclusion that development of generic total strain energy relationship as a 
measure of soil resistance against liquefaction by means of laboratory testing is feasible. 
If consistent sampling, sample handling, and testing methods and specifications are 
followed, the results are expected to be more consistent and vary within narrower limits. 
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Table 2.2 - Average Material Properties at the SRS Site 


PARAMETER/SOIL LAYER flocs 


cas 


SPT N-VALUE isis iets ss 
SHEAR WAVE VELOCITY, m/sec (fUs) 364 (1193) 381 (1251) 
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Test |.D.: MONT19 Controlled Parameter: Strain 
Relative Density (%): 60.6 Initial Effective Stress (kPa): 100 
Applied Shear Strain (%): 0.25 Frequency (Hz): 15 
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Figure 2.1 - Typical Time History Records of a Strain-Controlled Cyclic Triaxial Test 


Test I.D.: MONT19 Controlled Parameter: Strain 
Relative Density (%): 60.6 Initial Effective Stress (kPa): 100 
Applied Shear Strain (%): 0.25 Frequency (Hz): 15 
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Figure 2.2 - A Typical Plot of Shear Stress - Shear Strain Hysteresis Loops Developed 
During Cyclic Tnaxial Tests 
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Figure 2.3 - Grain Size Distribution Curve for Monterey No. 0 Sand (Arango, 1994) 
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Figure 2.4 - Strain Energy at Liquefaction Onset as a Function of Relative 
Density for Monterey No. 0 Sand - UCB Data 


2-20 


Energy at Liquefaction Onset, Eji, (J/m?*) 


5000 


4500 


4000 


3500 


3000 


2500 


2000 


1500 


1000 


500 


Fines Content = 2.0 percent 

2-Way Sinusoidal Loading 

Initial Effective Confining Pressure (o,0') = 100 kPa 
Monterey No. 0 Sand 


© Strain Controlled, Dr=41-42% 


@ Strain Controlled, Dr=61-62% 


O Stress Controlled, Dr=49-52% 


@ Stress Controlled, Dr=60-62% 


Frequency (Hz) 10 100 


Figure 2.5 - Strain Energy at Liquefaction Onset as a Function of Frequencies 
of Loading for Monterey No. 0 Sand - UCB Data 
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Figure 2.6 - Generalized Subsurface Soil Profile at SRS Site 
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Figure 2.7 - Typical Grain Size Distribution Curve for Tobacco Road Soil Material, SRS (Riemer and Seed, 1994) 
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Figure 2.8 - Strain Energy at Liquefaction Onset as a Function of Confining 
Pressure for SRS Soil Samples - UCB Data 
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Figure 2.10 - Strain Energy at Liquefaction Onset as a Function of Relative 
Density for Clean Sands - University of Colorado Data 
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Figure 2.11 - Strain Energy at Liquefaction Onset as a Function of Confining 
Pressure for Silty Sands - University of Colorado Data 
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Figure 2.12 - Grain Size Distribution Curves for the Northridge Site Soil Samples 
(Arango and Migues, 1996) 
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Figure 2.13 - Strain Energy at Liquefaction Onset as a Function of Relative 
Density for the Northridge Samples - UCB Data 
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Figure 2.15 - Strain Energy at Liquefaction Onset as a Function of Relative 
Density for Monterey No. 0 and Kasumigaura Sands - WSU 
Data 
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Figure 2.16 - Effect of Cyclic Loading Types on Strain Energy at 
Liquefaction Onset for Monterey No. 0 Sand - WSU Data 
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Figure 2.17 - Effect of Cyclic and Transient Loadings on Strain Energy at 
Liquefaction Onset for Monterey No. 0 Sand - WSU Data 
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Figure 2.18 - Grain Size Distribution Curves for Soil Samples from the Lower 
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San Fernando Dam and the Reid Bedford Sand (Figeroa et al., 1995) 
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Figure 2.19 - Comparison of Strain Energy at Liquefaction Onset for All 
Test Groups - Clean Sands 
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Figure 2.20 - Comparison of Strain Energy at Liquefaction Onset for All 
Test Groups - Silty Sands 
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CHAPTER 3 
WILDLIFE SITE: SOILS AND EARTHQUAKE DATA 


3.1 BACKGROUND 


The Parkfield and Wildlife seismologic instrumentation arrays are among the few 
liquefaction arrays in operation in the United States. Both arrays are maintained by the 
United States Geological Survey (USGS). The Wildlife instrument array was established 
in the fall of 1982 primarily to monitor pore pressure in sand deposits undergoing 
earthquake loading. The site was an obvious choice since this region is seismically active 
and the site has repeatedly liquefied during past earthquakes. The data recorded during 
one of the recent earthquakes (Superstition Hills earthquake of November 24, 1987) 
provide a unique opportunity to perform ground response analyses and to compare 
predicted ground motions with recorded motions. The results of these analyses in terms 
of the total strain energy can in turn, be used to compare the laboratory-derived strain 
energy with the energy calculated to cause liquefaction in the field. 


3.2 WILDLIFE SITE 


The Wildlife instrument array is located in the Imperial Valley, Southern California about 
36 km north of El] Centro (see Figure 3.1). This is one of the most seismically active 
regions in the United States. Recent major earthquakes occurred in this region in 1979, 
1981, and 1987. The array is on level ground and is located in the flood plain of the 
Alamo River, 20 m west of the river channel. Numerous articles, and several doctoral 
dissertations focus on the geological, geotechnical and seismological aspects of the site 
(Bennett et al., 1984; Bierschwale and Stokoe, 1984; Dobry et al., 1989; Etheredge et al., 
1987; Haag and Stokoe, 1985; Holzer et al., 1989; Magistrale et al., 1989; Matasovic and 
Vucetic, 1993; Vucetic and Thilakaratne, 1989; Youd et al., 1989; Wald et al., 1990; and 
Zeghal and Elgamal, 1994). Lateral spreading associated with liquefaction caused by the 
1987 Superstition Hills earthquake has been mapped by Holzer et al. (1989). The ground 
response analyses presented in this chapter was performed based on the data recorded 
during this event. 


3.33. STRATIGRAPHY AND SOIL PROPERTIES AT THE WILDLIFE SITE 


The subsurface stratigraphy at the Wildlife Site is shown in Figure 3.2 (Bennett et al., 
1984). Soil properties in terms of soil classification, index properties, penetration data 
(CPT and SPT), water content, and fines content for each subsurface soi] unit are shown 
in Figure 3.3. The subsurface units are defined as follows (Holzer et al., 1989; Bennett et 
al., 1984): 
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e Unit A is 2.3 m thick and consists of fine-grained flood plain deposits. The clay 
content in this layer is relatively high, about 16% to 65%. The water table is near the 
bottom of this unit, typically at a depth of 1.5 m. 


e Unit B1 is a loose, moderately sorted sandy silt about 1 m thick and highly 
susceptible to liquefaction (note the blow count of 3 in Figure 3.3 for this unit). The 
grain size distribution from the sand boil materials in the 1981 magnitude 5.9 
Westmoreland earthquake matched the gradation of the materials from this unit (Haag 
et al., 1985). 


e Unit B2 is a coarser deposit of silty fine sand to very fine sand about 3.5 m thick and 
is moderately dense. 


e Unit C consists of deposits of silty clay and clayey silt and is about 5 m thick. This 
unit and the layers below are very stiff and are not susceptible to liquefaction. 


Both seismic crosshole and the Spectral Analyses of Surface Waves (SASW) techniques 
have been used to measure the shear wave velocities at the site (Bierschwale and Stokoe, 
1984). The shear wave velocity profiles are shown in Figure 3.4. These profiles have 
been used in the ground response analyses performed for the current study as discussed 
below. 


Strain-dependent soil properties were obtained by Ladd (1982) and Turner and Stokoe 
(1982) from tests on soil samples taken at the Herber Road site, located about 40 km from 
the Wildlife Site. These soils are believed to consist of similar materials as those at the 
Wildlife Site. Shear modulus reduction curves and damping curves are shown in Figures 
3.5, 3.6, and 3.7. The values of G/Gmax and damping shown on the curves were selected 
at 10 cycles of loading to take into account the strain softening due to an earthquake 
loading of approximately 10 equivalent cycles (earthquake magnitude of about 6.6). 


3.4 INSTRUMENTATION AT THE WILDLIFE SITE 


The instrumentation at the Wildlife Site has been documented by Holzer et al. (1988) and 
is Shown in Figure 3.8. The instrumentation consists of: 


e Six pore pressure transducers installed in borings located along the circumferences of 
a 9.2 m diameter circle, with a horizontal spacing of 4.6 m. Two transducers are in 
Unit B1 (at 2.8 m and 3.0 m depths), three in Unit B2 (at 4 m, 5 m, and 6.3 m depths), 
and one in Unit D at 12.1 m. 


e Two triaxial accelerometers located in the center of the circle. The downhole 
accelerometer was placed at a depth of 7.5 m in Unit C, and a surface accelerometer 
was mounted on a concrete pad. 
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3.5 RECORDED EARTHQUAKE DATA 


On November 24, 1987, the Elmore Ranch earthquake with a surface wave magnitude of 
6.2 occurred at 0154 GMT. No liquefaction or significant pore pressure build-up was 
observed during this event. Approximately 11 hours later, a second earthquake, the 
Superstition Hills earthquake, occurred at 1315 GMT with a surface wave magnitude of 
6.6. During this event, massive liquefaction and sand boils followed by lateral spreading 
took place at the site. Instrumentation readings recorded during this event are described 
below. The epicentral distance and the duration of these two events are compared in 
Table 3.1. The location of the epicenters, 1987 surface rupture on the Superstition Hills 
fault, and their positions relative to the Wildlife Site are shown in Figure 3.9. 


Following the Superstition Hills earthquake, major ground surface cracks oriented in the 
east-west direction and displacements towards the Alamo river were mapped. Maximum 
displacements up to 200 mm (8 inches) were observed (Holzer et al., 1989). The records 
of all pressure transducers and accelerometers were obtained and processed. As indicated 
in Table 3.1, the 360 degree component of the records is in the north-south direction. The 
90 degree component is in the east-west direction. The acceleration, velocity, and 
displacement time histories of the record in the 360 degree direction from both the 
surface and the 7.5 m deep accelerometers are shown in Figure 3.10. The corresponding 
records for the 90 degree direction are shown in Figure 3.11. The vertical records are 
shown in Figure 3.12. 


Although the records indicate low frequency characteristics of the ground surface motion 
towards the tail of the record, identification of the initial liquefaction from the time 
history records is difficult. The acceleration response spectra of the records at 5% 
damping calculated from the records at the surface and at 7.5 m depth are compared in 
Figures 3.13 and 3.14, for horizontal and vertical components, respectively. Both the 
horizontal components depicted in Figure 3.13 show amplification of motion within the 
7.5 m soil thickness, with the main amplification taking place in the frequency range of 
0.5 Hz to 3 Hz. All records show that the earthquake waves were polarized and the 
intensity of the motion in the 360 degree direction was much larger than the intensity in 
the 90 degree direction. The vertical motion at the ground surface shown in Figure 3.14 
indicates a significant amplification (about 8 times) at about 20 Hz which cannot be 
explained. 


In a study of liquefaction at the Wildlife Site, Matasovic et al. (1993) analyzed the pore 
pressure records and determined the soil unit weight and the location of the water table 
are consistent with the pore pressure records. The pore pressure ratios obtained from the 
records of the four pressure transducers are shown in Figure 3.15. The records show that 
the P5 transducer (see Figure 3.8) reached the onset of liquefaction first, and that the 
onset of liquefaction of the other three transducers followed. The onset of liquefaction 
occurred at about 35 seconds after the start of the earthquake at transducer P5 and at 
about 90 seconds at the other three transducers. This may be due to the procedures 
followed in installation of the transducers and the densification of the soils at the 


3-3 


locations of the transducers. Holzer et al. (1988) reported that the transducers in their 
housings were pushed into the soil from the bottom of the bore hole and that the bore 
hole was filled with bentonite-cement grout, which could have disrupted hydraulic 
continuity and caused a delay in recording the maximum pore pressure. For transducers 
at the greater depths (P3 and P1), the delay in recording of the maximum pore pressure 
may also have been caused by the time required for migration of the express pore 
pressure from Unit B1, in which liquefaction occurred first, to the deeper depths in 

Unit B2. 
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Figure 3.1 - Map of California with Locations of Parkfield and Wildlife Liquefaction 
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Figure 3.2 - Wildlife Liquefaction Array - Geotechnical Properties (Bennett et al., 1984) 
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Figure 3.4 - Comparison of Shear Wave Velocity Profiles from Crosshole and SASW 
Tests at Wildlife Site (Bierschwale and Stokoe, 1984) 
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Figure 3.8 - Instrumentation of the Wildlife Site (Bennett et al., 1984) 
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Figure 3.9 - Location Map of the Wildlife Site and the Epicenters of the Elmore Ranch (M, 
= 6.2) and Superstition Hills (M, = 6.6) Earthquakes (Porcella et al., 1987) 
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Figure 3.10 - Earthquake Time Histories at 1315 GMT, November 24, 1987 at Wildlife 
Array - Horizontal Motions in 360° Direction 
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Figure 3.11 - Earthquake Time Histories at 1315 GMT, November 24, 1987 at Wildlife 
Array - Horizontal Motions in 90° Direction 
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Figure 3.12 - Earthquake Time Histories at 1315 GMT, November 24, 1987 at Wildlife 
Site - Vertical Motions 
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CHAPTER 4 


GROUND RESPONSE ANALYSES AND COMPARISON 
WITH THE LABORATORY DATA 


4.1 METHODS OF ANALYSES 


To evaluate the strain energy induced in the soil layers by the 1987 Superstition Hills 
earthquake, a series of ground response analyses were performed using the computer 
programs SHAKE (Schnabel et al., 1972), the modified program BDESRA (DESRA, Lee 
and Finn, 1978), and the direct computation of the stress and strain time histories from 
the motions recorded at the top and bottom of the liquefied layer. The modifications to 
the program DESRA were performed to incorporate the Martin-Dovidenkov soil model 
so that the user’s specified strain-dependent soil properties can be used. In the ground 
response analyses, the effects of variations of the shear wave velocity and the strain- 
dependent soil properties on the ground responses are evaluated. 


The computer programs SHAKE and DESRA are widely used for ground response 
analysis. The computer program SHAKE works in frequency domain and uses the 
equivalent linear method to handle soil nonlinearity. The computer program DESRA 
works in time domain and uses the hyperbolic stress-strain soil model to handle soil 
nonlinear effects. In each method, the shear stress and the shear strain response time 
histories were obtained at several depths. The time history of the strain energy at time t, 
E(t), was obtained from the following relationship: 


T 
E(T) = > Z(t). Ay (t) (4.1) 


where T(t) is the average shear stress from time t to t + At, and Ay(t) is the increment of 
the shear strain from time t to ttAt. The summation is computed for the total duration of 
the motion. Due to the one-dimensional formulation of the methodologies adopted for 
the ground response analyses, strain energy time histories at several elevations in both 
360 degrees and 90 degrees directions were first calculated separately. These two time 
histories at each elevation were subsequently combined at each time step to obtain the 
total strain energy time history at the respective elevation. All acceleration time histories 
used in the analytical study described in this chapter were from the November 24, 1987, 
1315 GMT Superstition Hills earthquake records. 


4.2 RESULTS OF THE GROUND RESPONSE ANALYSES USING THE 
COMPUTER PROGRAM SHAKE 


A representative soil profile was modeled for the Wildlife Site using the subsurface data, 
the shear wave velocity profiles, and the strain-dependent soil properties discussed in 
Chapter 3. The resulting SHAKE model is shown in Figure 4.1. The initial shear wave 
velocity profile used in this model is the average of the crosshole and the Spectral 
Analysis of Surface Waves (SASW) velocity profiles shown in Figure 3.4. 


The sensitivity of the results to the use of the crosshole, the SASW, or the average 
velocity profiles was first evaluated in three parallel SHAKE runs: the first with the 
crosshole velocity profile, the second with the SASW velocity profile, and the third with 
the average profile. The motion recorded at 7.5 m depth in the 360 degree direction was 
used as the input motion for the three runs. The responses at the ground surface from 
SHAKE output were then compared with the recorded surface motion. The comparisons 
are presented in Figures 4.2, 4.3, and 4.4, respectively, for the crosshole, SASW, and the 
average velocity profiles. The crosshole velocity profile significantly overpredicts the 
peak response. The agreement between the calculated and recorded spectral responses 
improve when the SASW velocity profile is used. However, generally a better agreement 
is obtained when the average velocity profile is used. 


Having determined the more representative velocity profile, the SHAKE deconvolution 
simulation was then studied. The soil profile used was the average velocity profile. The 
input motion used was the motion recorded at the ground surface. The response 
calculated at the depth of 7.5 m was compared with the motion recorded which is shown 
in Figure 4.5. This comparison shows a relatively good agreement in the frequency range 
of 1.5 Hz to 10 Hz. 


Third, the SHAKE simulation in the 90 degree direction was performed. The input 
motion used 7.5 m downhole recorded motion in the 90 degree direction. The soil 
properties used in this run were the strain-compatible properties obtained from the 
companion run using the 360 degree input motion. The calculated ground surface motion 
is compared with the recorded ground motions in Figure 4.6. The agreement, particularly 
at the peak response at 8 Hz, is very poor. However, it should be noted that as the pore 
pressure builds up, the ground response shifts to the low frequency range. The major 
contribution to the strain energy occurs at the time of the high pore pressure build-up at 
which time the soil stiffness is reduced significantly. Thus, the major attributes to the 
total strain energy are from the low frequency response of the soil layer. For this reason, 
good agreement at high frequencies, such as at 8 Hz for the 90 degree direction, may not 
be important for the purpose of the total strain energy computation. 


The results of analyses in terms of maximum acceleration, maximum shear stress, and 
maximum shear strain for each direction of shaking using the three velocity profiles are 
shown in Figures 4.7 through 4.12. Depending on the velocity profile and the earthquake 
component, the maximum estimated shear strain in the liquefiable sand layer ranges from 


0.25% to 0.6%. In general, the 90 degree component of the motion has a lower intensity 
and results in smaller stress and strain peak responses. 


Typical shear stress and shear strain time histories at the depth of 4.2 m obtained from the 
SHAKE analysis from the 360 degree input motion are shown in Figures 4.13 and 4.14, 
respectively. The corresponding hysteresis loops are shown in Figure 4.15. 


The time histories of the total strain energy at four depths obtained from the SHAKE 
responses using the 360 degree input motion are compared in Figure 4.16. The same 
results for the 90 degree input motion are shown in Figure 4.17. As shown in these two 
figures, the strain energy obtained from the 360 degree of input motion is approximately 
twice the energy using the 90 degree component of the input motion. A summation of the 
strain energy from both directions of shaking is shown in Figure 4.18. 


4.3 RESULTS OF THE GROUND RESPONSE ANALYSES USING THE 
COMPUTER PROGRAM BDESRA (MODIFIED DESRA) 


In order to evaluate the effects of modeling nonlinear soil behavior on ground responses, 
in addition to the program SHAKE, the computer program DESRA with total and 
effective stress options was also used. Due to the fully nonlinear formulation of this 
program, only the recorded motion at 7.5 m can be used as input motion. 


In the total stress analysis method, the strain-dependent soil properties shown in 
Chapter 3 were used in the modified DESRA (BDESRA) program. The input motion was 
specified at the 7.5 m depth and the response at the ground surface was compared with 
the recorded motion. The results are compared in Figures 4.19 and 4.20 using the 360 
degree and the 90 degree input motions, respectively. While the response comparisons 
show a better agreement in the frequency range of 2 Hz to 4 Hz in the 360 degree 
direction compared to the results obtained with the program SHAKE, the agreement in 
the 90 degree direction remains to be poor. The maximum shear stresses, and the 
maximum shear strains are shown in Figures 4.21 and 4.22, respectively. The time 
histories of the strain energy and the summation over both directions of shaking are 
shown in Figures 4.23, 4.24, and 4.25, respectively. 


The same analyses were repeated using the effective stress capability of BDESRA which 
is capable of computing the pore pressure build-up and its effect on the soil properties. 
The comparison of the responses at the ground surface from both the 360 degree and the 
90 degree input motion are shown in Figure 4.26 and 4.27, respectively. The results in 
the 360 degree direction under estimate the response in a wide frequency range. The 
agreement in the 90 degree direction is still very poor. The results in terms of the 
maximum shear stresses and the maximum shear strains are shown in Figures 4.28 and 
4.29, respectively. The maximum shear strain in the 360 degree direction reaches the 
peak value of about 1.7% (see Figure 4.29), well above the values computed by the other 
methods. The strain energy from each direction of shaking and the total energy are 
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shown in Figures 4.30, 4.31 and 4.32 , respectively. The strain energy computed by this 
method is about half the value computed by the other methods. In addition, the pore 
water pressure obtained from the effective stress DESRA analysis normalized to the 
effective overburden pressure at five depths have also been computed. The pore pressure 
ratios have been plotted in Figures 4.33 and 4.34 due to shaking from each direction. As 
shown in these figures, in the 360 degree direction, the pore pressure ratio reached the 
maximum value of one, while in the 90 degree direction it remains below the value one. 
In general, the results of the effective stress analyses are in poorer agreement with the 
recorded motion when compared to the total stress analyses results. This may be due to 
the fact that the strain-dependent soil properties used already includes the effect of the 
pore pressure build-up so that the use of the effective stress option double-counts the 
softening effect. 


4.44 RESULTS OF THE GROUND RESPONSE ANALYSES BASED ON 
DIRECT INTERPOLATION OF RECORDED MOTIONS 


In this approach, the recorded acceleration time histories at the top and bottom of the 
liquefied sand layer and the associated displacement time histories are used to estimate 
the average shear stresses and strains in the liquefied sand layer, assuming a linear 
distribution of the displacement and acceleration amplitudes over the thickness of the 
layer at each time step. The formulation of this approach is shown in Figure 4.35. It is, 
in essence, a linear interpolation of accelerations and displacement over the thickness of a 
layer based on recorded motions at the top and bottom of the layer. This approach is 
adopted from Zeghal and Elgamal (1994). As shown in Figure 3.8, the recording 
accelerometers are located at the surface and at the depth of 7.5 m below the liquefied 
layer. Assumptions for linear distribution of the stress and strain amplitudes may over 
estimate the actual amplitudes in the liquefied layer. Using this approach, the shear stress 
and the shear strain time histories are computed at several depths. The time histories 
corresponding to the depth of 5 m and the associated hysteresis loops due to the 360 
degree response motions are shown in Figures 4.36, 4.37, and 4.38, respectively. The 
same results for the 90 degree response motion are shown in Figures 4.39 through 4.41. 
The time histories of the strain energy computed using the ground motion for each 
direction of shaking and the total strain energy are shown in Figures 4.42 through 4.44. 
The results in terms of total energy are generally larger than the results obtained from the 
SHAKE and DESRA ground response analyses. This is mainly due to the fact that the 
accelerometer at 7.5 m is well below the liquefied unit B] (see Figure 3.2) which results 
in over estimation of shear strain in this approach. 
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4.5 | RESULTS OF THE GROUND RESPONSE ANALYSES USING THE 
COMPUTER PROGRAM SHAKE AND EPRI SOIL CURVES 


In a recent study, the Electric Power Research Institute (EPRI) published a set of generic 
soil modulus degradations and damping curves applicable to sand and clay soils under 
dry and submerged conditions (EPRI, 1993). The shear modulus and the damping curves 
used for the Wildlife ground response analysis are compared with the EPRI saturated 
sand curves in Figures 4.45 and 4.46, respectively. The Wildlife soil curves correspond 
to the upper bound of the shear modulus curves and the lower bound of the damping 
curves. The EPRI curves were also used in a series of SHAKE analyses using the 
recorded motion at the 7.5 m as the input motion. The response motions at ground 
surface in 360 degree and 90 degree directions are compared with recorded motions in 
Figures 4.47 and 4.48, respectively. No significant improvement in agreement with the 
recorded motions can be observed from the use of the EPRI soil curves. The shear stress 
and the shear strain and the associated hysteresis loops at a depth of 4.2 m are shown in 
Figures 4.49, 4.50, and 4.51, respectively. The strain energy for each direction of 
shaking and the total energy are shown in Figures 4.52 through 4.54. 


4.6 COMPARISON OF THE RESULTS OF GROUND RESPONSE 
ANALYSES 


The results of the ground response analyses in terms of total strain energy, obtained from 
SHAKE (using both the Wildlife soil curves and the EPRI soil modulus degradation and 
damping curves) and DESRA (total and effective analyses), are compared with the results 
derived from the recorded motions in Table 4.1 and Figure 4.55. The effective confining 
pressure shown in this figure is the mean effective stress at the respective depth of the 
response calculation. The results show a consistent trend except at a confining pressure 
of about 40 kPa to 50 kPa. This pressure corresponds to unit B2 (see Sections 3.3 and 
3.4), which may not have reached the onset of liquefaction during the earthquake. In fact, 
the pore pressure ratios reported by Matasovic and Vucetic (see Figure 3.15) show that 
the pressure transducers P] and P3 (at depths of 5 and 6.3 m) did not reach the pore 
pressure ratio of one, and that the large pore pressure ratios were developed after about 
50 seconds, indicating a redistribution of pore pressure after the shaking. For these 
reasons, the data below the depth of 5 m (16.4 ft) have been removed from the results and 
the new summary plot is depicted in Figure 4.56. As shown in this figure, the results of 
the ground response analyses obtained from all the different methods show a consistent 
trend and close agreement. In general, use of the EPRI soil curves in the SHAKE 
analyses results in larger total energy, mainly due to the differences in the damping 
curves (see Figure 4.46 ). In general, the results from the SHAKE analyses (using 
Wildlife soil curves) and the results obtained from the recorded motions are in better 
agreement. 
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4.7 COMPARISON OF THE RESULTS OF GROUND RESPONSE 
ANALYSES WITH THE LABORATORY DATA 


The results of the total strain energy for silty sands obtained from the available laboratory 
data are shown in Figure 2.16. As shown in this figure, the only set of data that covers 
the range of confining pressures less than 100 kPa are the summary results from Figueroa 
et al. Since the Wildlife liquefied layer is at a confining pressure lower less than 50 kPa, 
the relationship by Figueroa et al. (Equation 2.3) was used to extrapolate the data points 
to the low confining pressures. The extrapolated data points are plotted with the results 
from the ground response analyses in Figure 4.57. As shown in this figure, the agreement 
between the laboratory data and the results from the ground response analyses is good in 
spite of the fact that the laboratory data are not from the soil samples at the Wildlife Site. 
In addition, no correction factor has been applied to the laboratory data to correct for the 
field condition. Such factors should be evaluated once site-specific laboratory data are 
obtained. 
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Table 4.1 Summary of Strain Energy Computation 
from Ground Response Analyses 
(November 24, 1987, 1315 GMT Earthquake) 


Method of 


Average of SASW- 
and Crosshole Field 
Measurements 


Average of SASW 
and Crosshole Field 
Measurements 


BDESRA Average of SASW 
(Effective Stress) and Crosshole Field 
Measurements 


Direct Interpolation 
of Record Motions 


Soil Properties 


Analysis Shear Wave Velocity Strain Dependence | Ground Surface 
Profile Relationships (m) 


Depth from 


Channel Fill Sand 
(Ladd, 1982) 


EPRI (1993) 


Hyperbolic, Parameters 
adjusted to fit Channel 
Fill Sand (Ladd, 1982) 


Hyperbolic, Parameters 
adjusted to fit Channel 
Fill Sand (Ladd, 1982) 


Total 
Strain Energy 


(J/m*) 


Note: The motions recorded at 7.5 m depth were used as input motion in ground response analyses. 
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Figure 4.1 - Wildlife Site Soil Profile Used in Response Analyses Based on Average 
of SASW and Crosshole Shear Wave Velocity Measurements 
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Figure 4.7 - Maximum Acceleration Distribution over Depth from SHAKE Results - 
Wildlife Site, NOvember 24, 1987, 1315 GMT Earthquake in 360° Direction 
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Figure 4.8 - Maximum Acceleration Distribution over Depth from SHAKE Results - 


Soil properties were based on the 
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shear wave velocity measurements. 
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using the stronger NS (360 degrees) 
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iteration. 
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Wildlife Site, November 24, 1987, 1315 GMT Earthquake in 90° Direction 
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Figure 4.9 - Maximum Shear Stress Distribution over Depth from SHAKE Results - 
Wildlife Site, November 24, 1987, 1315 GMT Earthquake in 360° Direction 
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Figure 4.10 - Maximum Shear Stress Distribution over Depth from SHAKE Results - 
Wildlife Site, November 24, 1987, 1315 GMT Earthquake in 90° Direction 
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Figure 4.11 - Maximum Shear Strain Distribution over Depth from SHAKE Results - 
Wildlife Site, November 24, 1987, 1315 GMT Earthquake in 360° Direction 
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Figure 4.12 - Maximum Shear Strain Distribution over Depth from SHAKE Results - 
Wildlife Site, November 24, 1987, 1315 GMT Earthquake in 90° Direction 
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Figure 4.15 - Shear Stress - Shear Strain Hysteresis Loop at Depth 4.21 m 
(13.8 ft.) Calculated from SHAKE Output - Wildlife Site, 
November 24, 1987, 1315 GMT Earthquake in 360° Direction. 
Soil Properties are Based on the Average of SASW and 
Crosshole Shear Wave Velocity Measurements 
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Figure 4.21 - Maximum Shear Stress Distribution over Depth from BDESRA Results of 
Total Stress Analyses - Wildlife Site, November 24, 1987, 1315 GMT 
Earthquake 
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Figure 4.22 - Maximum Shear Strain Distribution over Depth from BDESRA Results of 
Total Stress Analyses - Wildlife Site, November 24, 1987, 1315 GMT 
Earthquake 
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Figure 4.28 - Maximum Shear Stress Distribution over Depth from BDESRA Results of 
Effective Stress Analyses - Wildlife Site, November 24, 1987, 1315 GMT 
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Figure 4.29 - Maximum Shear Strain Distribution over Depth from BDESRA Results of 
Effective Stress Analyses - Wildlife Site, November 24, 1987, 1315 GMT 
Earthquake 
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Figure 4.38 - Shear Stress - Shear Strain Hysteresis Loop at Depth of 5.06 m (16.6 ft.) 
- Based on Direction Interpolation of Recorded Motions, Wildlife Site, 
November 24, 1987, 1315 GMT Earthquake in 360° Direction. 
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Figure 4.41 - Shear Stress - Shear Strain Hysteresis Loop at Depth of 5.06 m (16.6 ft.) 
- Based on Direct Interpolation of Recorded Motions, Wildlife Site, 
November 24, 1987, 1315 GMT Earthquake in 90° Direction 
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Figure 4.51 - Shear Stress - Shear Strain Hysteresis Loop at Depth of 4.21 m (13.8 ft.). 
Calculated from SHAKE Output - Wildlife SIte, November 24, 1987, 
1315 GMT Earthquake in 360° Direction. Soil Properties Are Based on 
the Average of SASW and Crosshole Shear Wave Velocity Measurements. 
The EPRI (1993) G/Gmax and Damping Curves Are Used in the Calculation. 
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CHAPTER 5 
SUMMARY AND RECOMMENDATION 


A feasibility study was conducted to determine the applicability of the strain energy 
method for liquefaction potential evaluation. This study consisted of collection and 
synthesis of available laboratory test data and evaluation of the strain energy as a measure 
of soil resistance to liquefaction. The data obtained by several groups of researchers 
confirm the observation that the strain energy (in normalized form) follows the time 
history of the pore pressure ratio increase up to the onset of liquefaction. The data also 
show that by using the relative density for clean sands and the confining pressure for silty 
sands, development of the generic strain energy relationships for liquefaction resistance is 
feasible. However, other factors such as the shape and size of the sand particles, the type 
and the amount of fines content, and the effect of the rate and type of loading on the 
strain energy need to be investigated and quantified. 


The feasibility study also successfully predicts the strain energy required for liquefaction 
to take place in the field based on the recorded data at the Wildlife Site. The ground 
response analyses performed for this site, using the conventional methods of the ground 
response analyses, resulted in consistent results, with a scatter typically expected in such 
analyses. However, all methods used in the ground response analyses were one- 
dimensional, and the summation of the strain energy due to shaking in each direction is at 
best an approximation. Nevertheless, the comparison of the laboratory data with the 
results of ground response analyses is encouraging. Such favorable comparison suggests 
also that the goal of developing the strain energy approach for generic application to the 
evaluation of the soil liquefaction potential should be pursued in future phases of the 
project. 


Based on the results of the feasibility study, it is recommended to: 


e Perform laboratory testing, preferably on soil samples from the Wildlife Site 
at confining pressures in the range of 26 kPa to 50 kPa to cover the low end of 
the test results. 


e Perform ground response analyses using a nonlinear computer program with 
multi-directional shaking capability to validate or modify the summation of 
the total strain energy calculated by the methods limited to uni-directional 
shaking capability as used in the current study. 


e Apply the energy based method to liquefaction/no liquefaction sites at 
Northridge, Kobe in Japan where recorded motions are available, and at the 
Wildlife Site for the Elmore Ranch earthquake during which no liquefaction 
occurrence was observed. 
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The energy method provides a sound approach to the evaluation of liquefaction potential. 
It is based on the basic principles of the “stress” and “strain” approaches in which the soil 
capacity to resist liquefaction is measured in terms of the strain energy. It can also handle 
the special characteristics of the ground motion, such as the near-field effects. The 
method also has the capability of providing an estimate of the amount of pore pressure 
build-up before the onset of liquefaction. In this approach, the laboratory data can also 
be corroborated as more field data from the sites that have or have not liquefied become 
available. 
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LABORATORY TESTS ON MONTEREY NO. 0 SAND, PERFORMED AT THE 
UNIVERSITY OF CALIFORNIA, BERKELEY 
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Test I.D.: MONT10 Controlled Parameter: Stress 
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Test I.D.: MONT11 Controlled Parameter: Stress 
Relative Density (%): 60.6 Initial Effective Stress (kPa): 100 
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Test I.D.: MONT12 Controlled Parameter: Stress — 
Relative Density (%): 61 initial Effective Stress (kPa): 100 
Applied Stress Ratio: 0.23 Frequency (Hz): 1 
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Test I.D.: MONT14 Controlled Parameter: Stress 
Relative Density (%): 60.2 Initial Effective Stress (kPa): 100 
Applied Stress Ratio: 0.3 Frequency (Hz): 20 
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Test I.D.: MONT17 Controlled Parameter: Strain. 
Relative Density (%): 61 Initial Effective Stress (kPa): 100 
Applied Shear Strain (%): 0.25 Frequency (Hz): 1 
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Test I.D.: MONT18 Controlled Parameter: Strain 
Relative Density (%): 61 Initial Effective Stress (kPa): 100 
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Test I.D.: MONT18 Controlled Parameter: Strain 
Relative Density (%): 61 Initial Effective Stress (kPa): 100 
Applied Shear Strain (%): 0.25 Frequency (Hz): 10 
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Test I.D.: MONT19 Controlled Parameter: Strain. 
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Test I.D.: MONT20 Controlled Parameter: Strain — 
Relative Density (%): 40.9 Initial Effective Stress (kPa): 100 
Applied Shear Strain (%): 0.25 Frequency (Hz): 1 


Shear Stress, + (kPa) 


time, t (sec) 


Shear Strain, y (%) 


time, t (sec) 


Energy per Volume, E (J/m’) 


0 2 4 6 8 10 12 
time, t (sec) 


6 
time, t (sec) 


d:\nit\itp\cy_filesWMONT20.XLS 8/15/96 11:47 AM 
A-21 


Test I.D.: MONT20 Controlled Parameter: Strain 
Relative Density (%): 40.9 Initial Effective Stress (kPa): 100 
Applied Shear Strain (%): 0.25 Frequency (Hz): 1 
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Test I.D.: MONT22 Controlled Parameter: Strain | 
Relative Density (%): 42.3 Initial Effective Stress (kPa): 100 
Applied Shear Strain: 0.08% Frequency (Hz): 1 
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Test I.D.: MONT24 Controlled Parameter: Stress 
Relative Density (%): 51.8 Initial Effective Stress (kPa): 100 
Applied Stress Ratio: 0.31 Frequency (Hz): 1 
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Test I.D.: MONT25 Controlled Parameter: Stress 
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Test I.D.: MONT26 Controlied Parameter: Stress 
Relative Density (%): 49.9 Initial Effective Stress (kPa): 100 
Applied Stress Ratio: 0.35 Frequency (Hz): 1 
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Test I.D.: MONT26 Controlled Parameter: Stress 
Relative Density (%): 49.9 Initial Effective Stress (kPa): 100 
Applied Stress Ratio: OFS Frequency (Hz): 1 


Shear Stress vs. Shear Strain 


t 
pay 
(eo) 


Shear Stress, + (kPa) 


Shear Strain, y (%) 


d:\nit\itp\cy_files\MONT26.XLS 8/21/96 8:58 AM 
Aq 32 


Test I.D.: MONT30 Controlled Parameter: Strain’ 
Relative Density (%) 41.9 Initial Effective Stress (kPa): 100 
Applied Shear Strain (%): 0.26 Frequency (Hz): 10 


Shear Stress, + (kPa) 


time, t (sec 


tel 5 ] ] j 
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
1 
“02 


time, t (sec) 


Shear Strain, y (%) 


8 & 


N 
rs) 


Energy per Volume, E (J/m’) 
§ 
rs) 


oO 


0) 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
time, t (sec) 


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
time, t (sec) 


d:\nit\itp\cy_files\MONT30.XLS 8/21/96 9:08 AM 
A-33 


Test I.D.: 
Relative Density (%) 
Applied Shear Strain (%): 


Shear Stress, + (kPa) 


0.4 -0.3 


d:\nit\itp\cy_files|MONT30.XLS 


MONT30 Controlled Parameter: 


41.9 Initial Effective Stress (kPa): 


0.26 Frequency (Hz): 


Shear Stress vs. Shear Strain 


3 
Shear Strain, y (%) 


A-34 


Strain — 
100 
10 


0.3 


8/21/96 9:08 AM 


MONT33 


— ul 
«| i i 
Bel ATTA 


= Ce ee ee 


: | —_—_—_________., 


———————S 
——————— 
————————} 
——————— 
ee 
———_} 
——J 
es 


a 


—~———— —— 
S——S 


—————— 


—— 
——_— 
= : 
| ——= 
£ 


Test |.D.: MONT33 Controlled Parameter: Strain 
Relative Density (%) 40.5 Initial Effective Stress (kPa): 100 


Applied Shear Strain (%) 0.08 Frequency (Hz): 10 
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Test I.D.: MONT35 Controlled Parameter: Stress 
Relative Density (%) 61.8 Initial Effective Stress (kPa): 100 
Applied Stress Ratio: 0.35 Frequency (Hz): 0.1 
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Test I.D.: MONT35 Controlled Parameter: Stress, 
Relative Density (%): 61.8 Initial Effective Stress (kPa): 100 
Applied Stress Ratio: 0:35 Frequency (Hz): 0.1 
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Test I.D.: MONT37 Controlled Parameter: Strain 
Relative Density (%): 61.6 Initial Effective Stress (kPa): 100 
Applied Shear Strain (%) On Frequency (Hz): 1 
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Test I.D.: MONT37 Controlled Parameter: Strain 
Relative Density (%): 61.6 Initial Effective Stress (kPa): 100 
Applied Shear Strain (%) 0.17 Frequency (Hz): 1 
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Test I.D.: MONT38 Controlled Parameter: Strain. 
Relative Density (%): 61.8 Initial Effective Stress (kPa): 100 
Applied Shear Strain (%): 0.25 Frequency (Hz): 10 
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Test I.D.: MONT38 Controlled Parameter: Strain 
Relative Density (%): 61.8 Initial Effective Stress (kPa): 100 
Applied Shear Strain (%): 0.25 Frequency (Hz): 10 
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Test I.D.: B23P3BCY Controlled Parameter: Stress. 
Fines Content (%): 16.6 Initial Effective Stress (kPa): 200 
Dry Density (kN/m*): 16.4 


ON 


time, t (sec) 


——— 
—— 


——— or 
—_—_—_—— 


Shear Stress, + (kPa) 


Shear Strain, y (%) 


-15 

time, t (sec) 

—~ 80000 

= 

= 

uw 60000 

a 

= 

= 

% 40000 

> 

g 

> 20000 

om 

a 

i 

_ 04 

0 20 40 60 80 
time, t (sec) 
40 
time, t (sec) 
d:\nit\itp\cy_files\B23P3BCY.XLS 8/21/96 7:48 AM 


B~9 


Stress 
200 


WW 
WE 
Wy 


sini 
A 
I | 

| Wa WHA A 


| 


SS 


| 


i i 
ih I 
mi 


— 


| 
/ 
y! 
ih 
! | 
ity 
| 


Controlled Parameter: 
Initial Effective Stress (kPa): 


B23P3BCY 


16.6 
16.4 


Content (%): 


Test I.D.: 
Fines 


Dry Density (kKN/m*): 


Shear Strain 


Shear Stress vs. 


d:\nit\itp\cy_files\B23P3BCY.XLS 


Test I.D.: B23P3MCY Controlled Parameter: Stress. 
Fines Content (%): 18.5 Initial Effective Stress (kPa): 200 
Dry Density (kKN/m*): 15.22 


Shear Stress, < (kPa) 


Shear Strain, y (%) 


time, t (sec) 


20000 


10000 


Energy per Volume, E (J/m’) 


time, t (sec) 


Ru E/Eliq 


time, t (sec) 


d:\nit\itp\cy_files\B23P3MCY.XLS 8/21/96 7:52 AM 
B=11 


B23P3MCY 


——— 


WwW 
oO 
SS SSS 
ee 


h 
nN N 
a 
——SS=_="= Soe 
Ss oceans : 
S a 


I 


— 


262. SSS = 
= <= 
——— 
nF da 


—s 


\ 
WS 
Sy 
1S Se 
SSS 
—= 
Se == 


— 
== Stes 
nes 
—— 


—SS 


SS 


Se,’ 


—s 


\N 
rs 
a an Sn 
Pk 


_enDe = eeaeee 
SSG ES ——— 
ss , 
Sn 
se 


FSS 


MY 
| 

| 
OF 


Test I.D.: B23P3TCY Controlled Parameter: Stress 


Fines Content (%): 20.5 _ Initial Effective Stress (kPa): 200 
Dry Density (kN/m‘*): 16.63 
50 
40 | 
zx | 
E 2 | | 
i 
® fe) | | 
® -10 0 5 10 15 20 25 
S 20 | 
® -30 | | | 
-40 \ f 


time, t (sec) 


Shear Strain, y (%) 


time, t (sec) 


: 


: 


20000 


Energy per Volume, E (J/m’) 


0 cS) 10 time, t (sec) 15 20 25 


NY ae A Sli. 


d:\nit\itp\cy_files\B23P3TCY.XLS 8/21/96 7:57 AM 
B-13 


d:\nit\itp\cy_files\B23P3TCY.XLS 


wil 
" 


vn 


a 


a 
—— 
_————— 


ial 


5 UU NE nl 
ULES esl ii TAA i | ii i | y 


time, t (sec) 


Test I.D.: 
Fines Content (%): 
Dry Density (kN/m’): 


Shear Stress, +.(kPa) 


d:\nit\itp\cy_files\B12PSBCY.XLS 


B12P5BCY Controlled Parameter: 


27 Initial Effective Stress (kPa): 


15.82 


Shear Stress vs. Shear Strain 


Shear Strain, y (%) 


B-16 


Stress" 
300 


8/20/96 4:34 PM 


d:\nit\itp\cy_files\B12PSMCY.XLS 


rest 1.D:: Bi2P5MCY Controlled Parameter: Stress 
Fines Content (%): 26.8 Initial Effective Stress (kPa): 300 
Dry Density (kN/m*): 16.78 


EET 
aA TT 


time (sec) 


; ‘bani ta 
: : VV Tol YY NANI titi 


-10 


ee 


SS SSS 
———————— 


Shear Stress, + (kPa) 
Sees 68s 3 8 
(2) a = 


= 8 


= = 
(=) (oa) 


Sete 
gee 
———_ 
=> 
See) 
ea 
—— 
pie SSS 


See 
= 


=== 8 


time, t (sec) 


350000 


200000 


0 20 time, t (sec) a0 om 


0 10 20 30 40 50 60 70 
time, t (sec) 


B-17 


8/20/96 4:42 PM 


Stress. 
300 


Controlled Parameter: 


Bi2P5MCY 


26.8 


Test I.D.: 


w 
eas 


= 
SoS 


Initial Effective Stress (kPa) 


“ 


Xe 
IW 


Shear Stress vs. Shear Strain 


16.78 


Fines Content (%) 
Dry Density (kN/m*) 


(2qy) > ‘ssaajg seays 


Shear Strain, y (%) 


8/20/96 4:42 PM 


d:\nit\itp\cy_files\B12P5MCY.XLS 


B-18 


Test I.D.: B12P5TCY Controlled Parameter: Stress — 


Fines Content (%): 22.3 Initial Effective Stress (kPa): 300 
Dry Density (kN/m*): 18.02 
80 4 
60 
< 40 | 
=< 
© 20 | 
a 
@ 0 
® 9 0 20 40 
o 
= -40 
w” 
-60 
ie yy VV 
time (sec) 


Shear Strain, y (%) 


time, t (sec) 


150000 


100000 


Energy per Volume, E (J/m’) 


0 20 time, t (sec) 40 


Wires salle | 


0 20 time, t (sec) 40 


d:\nit\itp\cy_files\B12PSTCY.XLS 8/20/96 4:45 PM 
B-19 


Test I.D.: B12P5TCY Controlled Parameter: S) 
Fines Content (%): 22:3 Initial Effective Stress (kPa): 300 


Dry Density (kN/m*): 18.02 


Shear Stress vs. Shear Strain 


7 


| 


MTAIVMAA NINA 
7 


| 


———— 


Ss CO 

SN 
SS ——————————— 
_——— 


f 
j 
( 


—— 


| 


———————— 


| 


SSS SSS 
SSS 


SS 
— 
—S 


Wa 


My 


a 


a 
2s 
= 3 
—— 


——— 


Hn 


i 


| 


, 


Wi 


HH 


d:\nit\itp\cy_files\B12PSTCY.XLS 


Test I.D.: B12P7BCY Controlled Parameter: Stress © 
Fines Content (%): 15 ch Initial Effective Stress (kPa): x Ff 
Dry Density (kN/m’): 15.91 


Shear Stress, + (kPa) 


Shear Strain, y (%) 


20 0 10 20 30 


time (sec) 


ey) 


time, t (sec) 


Energy per Volume, E (J/m’) 
Beg: 
so: 5 


0 10 time, t (sec) 20 30 


15.0 


10.0 


5.0 


0 10 time, t (sec) 20 30 


d:\nit\itp\cy_files\B12P7BCY.XLS Bo 8/20/96 4:50 PM 


Test I.D.: B1i2P7MCY Controlled Parameter: Stress 


Fines Content (%): Ly. Initial Effective Stress (kPa): CWS) 
Dry Density (kN/m’): 16.25 
: | 

_ 4 

iJ 

= 20 

_ 

é 

7 

oS 

= 

” 


time (sec) 


Shear Strain, y (%) 
Oo 


5 
-10 
-15 : 
time, t (sec) 

150000 
g 
= 
wi 

@ 100000 
E 
= 
o 
> 

& 50000 
> 
= 
) 
i= 
Ww 

0 
0 10 20 time, t (sec) 30 40 50 


0 10 20 ~='time,t (sec) 30 40 50 


d:\nit\itp\cy_files\B12P7MCY.XLS - 8/20/96 4:54 PM 
B- 


Bi2P7MCY 
17 


| vi i) 
: y Bd 


A i ra | 
Wh 


i | | 


! 4 | 


ie 


SS 


Test I.D.: B12P7TCY Controlied Parameter: Stress 
Fines Content (%): Adi, Initial Effective Stress (kPa): oN he) 
Dry Density (kN/m’): > 16.24 


80 
60 
| 

40 | 
20 

: Baa | 
2.0 5 10 15 20 
-40 
60 
-80 


Shear Stress, + (kPa) 


time, t (sec) 
10 
<5 
nd 
£ 
= 0 | 
_ 0 5 10 15 20 
® 
= 5 
-10 
time, t (sec) 
=~ 40000 
E 
= 
wi 30000 
o 
£ 
= 
%o 20000 
> 
= 
>, 10000 
= 
+) 
on 
ld 04 
0) S 10 15 20 
time, t (sec) 


10.0 


5.0 


0 5 10 15 20 
time, t (sec) 


d:\nit\itp\cy_files\B12P7TCY.XLS 8/21/96 9:48 AM 
B-25 


Test I.D.: 


Fines Content (%): 
Dry Density (kN/m’): 


<-(kPa) 


Shear Stress, 


d:\nit\itp\cy_files\B12P7TCY.XLS 


B1i2P7TCY 
15.7 
16.24 


Shear Stress vs. Shear Strain 


Controlled Parameter: 


Initial Effective Stress (kPa): 


Shear Strain, y (%) 


Stress 
S15 


8/21/96 9:48 AM 
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Test I.D.: B2P6MCY Controlled Parameter: Stress 
Fines Content (%): 20.1 Initial Effective Stress (kPa): 743 
Dry Density (kN/m*): 16.5 


Shear Stress vs. Shear Strain 


200 


Shear Stress, < (kPa) 


Shear Strain, y (%) 


d:\nit\itp\cy_files\B2P6MCY.XLS 8/21/96 8:33 AM 
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APPENDIX C 


LABORATORY TESTS ON SOIL SAMPLES, PERFORMED AT THE 
UNIVERSITY OF COLORADO 


Aem-z jeplosnuls 
Aem-z jeplosnuls 

eBM-z JBEPIOSNUIS 
Aem-z jepiosnuls 

eM-Z jePlOSNulS 
Aen-z |e piosnuls 
Aem-z jepilosnuls 
Aem-Z jeplosnuis 


eM-z |BEPlosnuls 


adeys peo} 


[Tear 7 [ Tonues 


‘payoeyy Ue Sg] ISeL JO aduanbas ay} Bulmojjoy syjnsey \so1 


v94 
vos 
9rd 
9b 
1 


€Losaon 


6040N 
ZOA0N 


OPRIOJOD JO AjIsIOAIUA 1e paULIOjiog spueg AjIg pue UI]_D UO BIG ISI] [BUOISIOT, NDAD oy} Jo ArBUILUNS - [°- 91qeL, 


C-2 


Test I.D.: UOFCS Controiled Parameter: Stress 
Fines Content (%): Clean Sand Initial Effective Stress (kPa): 200 
Relative Density (%): 32.6 


Shear Stress, + (kPa) 


time, t (sec) 


Shear Strain, 7 (%) 
ro) 
8 
& 
3 
3 
8 


time, t (sec) 


Energy per Volume, E (J/m’) 


0 20 40 time, t (sec) 60 80 100 


0 20 40 time, t (sec) 60 80 100 


d:\nit\itp\cy_files\UOFC5.XLS 8/19/96 4:18 PM 
C-3 


Test I.D.: UOFCS5 Controlled Parameter: Stress’ 
Fines Content (%): Clean Sand Initial Effective Stress (kPa): 200 
Relative Density (%): 32.6 
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Test I.D.: UOFC7 Controlled Parameter: Stress © 
Fines Content (%): Clean Sand Initial Effective Stress (kPa): 205 
Relative Density (%): 41 
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Test I.D.: UOFC7 Controlled Parameter: Stress 
Fines Content (%): Clean Sand Initial Effective Stress (kPa): 205 
Relative Density (%): 41 
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Test I.D.: UOFC9 Controlled Parameter: Stress 
Fines Content (%): Clean Sand Initial Effective Stress (kPa): 304 
Relative Density (%): 42 
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Test I.D.: UOFC9 Controlled Parameter: Stress © 
Fines Content (%): Clean Sand Initial Effective Stress (kPa): 304 
Relative Density (%): 42 
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Test I.D.: UOFC13 Controlled Parameter: Stress 
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Test I.D.: UOFC13 Controlled Parameter: Stress 
Fines Content (%): 20 Initial Effective Stress (kPa): 300 
Dry Density (kN/m*) 15.26 
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Test I.D.: UOFC14 Controlled Parameter: Stress’ 
Fines Content (%): 20 Initial Effective Stress (kPa): 200 
Dry Density (kN/m*) 15.18 
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Test I.D.: UOFC14 Controlled Parameter: Stress’ 
Fines Content (%): 20 Initial Effective Stress (kPa): 200 
Dry Density (kN/m’*) 15.18 
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Test I.D.: UOFC15 Controlled Parameter: Stress 
Fines Content (%): 20 Initial Effective Stress (kPa): 200 
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Test I.D.: UOFC15 Controlled Parameter: Stress 
Fines Content (%): 20 Initial Effective Stress (kPa): 200 
Dry Density (kN/m*) 15.25 
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Test I.D.: UOFC17 Controlled Parameter: Stress © 
Fines Content (%): 45 Initial Effective Stress (kPa): 203.4 
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Test I.D.: UOFC17 Controlled Parameter: Stress 
Fines Content (%): 45 Initial Effective Stress (kPa): 203.4 
Dry Density (kN/m*) 16.18 
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Test I.D.: UOFC18 Controlled Parameter: Stress’ 
Fines Content (%): 45 Initial Effective Stress (kPa): 190 
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Test I.D.: UOFC18 Controlled Parameter: Stress 
Fines Content (%): 45 Initial Effective Stress (kPa): 190 
Dry Density (kN/m’) 16.18 
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APPENDIX D 


LABORATORY DATA ON SOIL SAMPLES FROM THE NORTHRIDGE SITE, 
PERFORMED AT THE UNIVERSITY OF CALIFORNIA, BERKELEY 
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Test I.D.: BTC3CY1 Controlled Parameter: Stress 
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Test I.D.: BTC3CY1 Controlled Parameter: ‘Stress 
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Test I.D.: BTC3CY2 Controlled Parameter: Stress 
Fines Content (%): 5 Initial Effective Stress (kPa): 100 
Dry Density (kN/m’) 16.12 
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Test I.D.: BTC4CY1 Controlled Parameter: Stress 
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Test I.D.: BTC4CY1 Controlled Parameter: Stress 
Fines Content (%): 5 initial Effective Stress (kPa): 100 
Dry Density (kN/m*) 16.35 
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Test I.D.: BTC4CY2 Controlled Parameter: Stress. 
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Test I.D.: BTC6CY1 Controlled Parameter: Stress 
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Test I.D.: BTC6CY1 Controlled Parameter: Stress © 
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